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Finite energy QCD sum rules involving both inverse and positive moment integration kernels 
are employed to determine the bottom quark mass. The result obtained in the MS scheme at a 
reference scale of lOGeV is mi,(10GeV) = 3623(9) MeV. This value translates into a scale invariant 
mass rni,{mi,) = 4171(9) MeV. This result has the lowest total uncertainty of any method, and is 
less sensitive to a number of systematic uncertainties that affect other QCD sum rule determinations. 
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I. INTRODUCTION 



With the availability of new cross section data on e+e~ 
annihilation into hadrons from the BABAR collaboration 
[l[ , the bottom quark mass was determined recently with 
unprecedented precision using inverse moment QCD sum 
rules • The result in the MS scheme at a reference scale 
of 10 GeV is 



TOfc(10 GeV) = 3610(16) MeV . 



(1) 



However, as was subsequently pointed out [3], this result 
relies on the assumption that PQCD is already valid at 
the end point of the BABAR data, i.e. a/s = 11.21 GeV, 
where s is the squared energy. This assumption might be 
questionable, as the prediction of PQCD for the i?-ratio 
does not agree with the experimentally measured value at 
this point. This QCD sum rule result was also shown to 
depend significantly on this assumption. Hence, further 
reductions in the error of the bottom-quark mass using 
QCD sum rules will depend on the ability to control this 
systematic uncertainty. One way of achieving this would 
be for a new experiment to extend the BABAR measure- 
ment into a region where PQCD is unquestionably valid. 
In this paper we follow another approach based entirely 
on theory. We use a finite energy QCD sum rule with 
integration kernels involving both inverse and positive 
powers of the energy, as employed recently to determine 
the charm-quark mass We also exploit the freedom 
offered by Cauchy's theorem to reduce the dependence 
of the quark mass on the above systematic uncertainty. 
This is achieved by using integration kernels that re- 
duce the contributions in the region y/s ~ 11.21 GcV to 
^/so, where there is no data and the onset of PQCD at 



s — So has to be assumed. As a benefit, this procedure 
reduces also the continuum contribution relative to the 
well known T narrow resonances. 



II. THEORETICAL BACKGROUND 



We consider the vector current correlator 



n^.(g') = * I d^xe'i^O\T{V^ix)V,M)\0) 



(2) 



where V^{x) = b{x)^^b(x), and h{x) is the bottom-quark 
field. Cauchy's residue theorem in the complex s-plane 
{—q^ = Q"^ = s) implies that 

1 If 

/ p{s)—lmll{s)ds — (b p{s)Il{s)ds 

Jo I" 27ri Jc(\so\) 

+ Res[H(s)p(s),s = 0], (3) 
where p{s) is an arbitrary Laurent polynomial, and 

1 



ImH(s) = —Rbis) 
iZn 



(4) 



with Rb{s) the standard i?-ratio for bottom production. 
The power series expansion of n(s) for large and space- 
like s can be calculated in PQCD, and has the form 

n(^)|p<,cB-?E(^)"n'"'W, (5) 

n=0 
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where e;, = 2/3 is the bottom-quark electric charge, and 

n(«)(,).^(^)'n("). (6) 

Here rub = rni,{fj.) is the quark mass in the MS scheme, 
at the renormaUzation scale /i. The order Ola^iml/ sY] 
results for i = 1, • • • , 6 have been calculated in ^Sj, with 
new results up to OlaKml / s)^^] obtained recently Q. 

At order 0[al], n[,^^ and nf^ are known 0], and the 

logarithmic terms in XIj"^-* may be found in Q. The con- 

(3) 

stant term in '^ot known exactly, but has been 

estimated using Pade approximants @, and the Mellin- 
Barnes transform fTo'|. At order 0[af] the exact logarith- 
mic terms in nQ"*-* and 11 were determined in [Tl|-[T^. 
whilst the constant terms are not yet known. Given that 
these constant terms will contribute to sum rules with 
kernels containing powers and s*^, respectively, for 
consistency we shall not include any five-loop order ex- 
pressions. However, we find that if all known five-loop or- 
der terms are taken into account, the mass of the bottom- 
quark only changes by roughly 0.03%, which is about a 
tenth of the accuracy of this determination. 
The Taylor series expansion of H(s) about s = is usu- 
ally cast in the form 

n(^)|pQCD=|iE^»-"' (7) 

n>0 

where z = s/ (4to^). The coefficients Cn can be expanded 
in powers of as (fJ-) as 




where = ln(m^//.t^). Up to O ja^), the coefficients up 
to n = 30 of Cn are known [l3|-[3|- There is also a sub- 
leading contribution of order 0{a^ (jric/frib)^) [3 affect- 
ing the coefficient C„ in Eq.®, as well as QED correc- 
tions. The former contributes around — l.OMeV, and the 
latter roughly -2.0 Me V to the resuh for mb ( 10 GeV). 
Finally, there is a non-perturbative contribution to H(s) 
from the gluon condensate, but it has been found to be 
completely negligible [11]. We fully agree, and thus con- 
firm this result. For the strong running coupling we use 
the Particle Data Group ^ value as{mz) = 0.1184(7), 
which corresponds to Q;^(10GeV) = 0.1792(16). 



0.6 




0.1 - 



10.6 10.7 10.8 10.9 11.0 11.1 11.2 11.3 

•/7 (GeV) 

FIG. 1: The corrected BABAR data 1] and the PQCD pre- 
diction (solid black line) obtained using Rhad ^20]. 



III. EXPERIMENTAL INPUT 



In order to evaluate the left-hand side of Eq. ^ one needs 
to use experimental input. First, there are the four nar- 
row T-resonanccs, and we calculate their contribution to 
Eq.Q using the zero- width approximation 

e.^^9W^^^^_ (9) 

where i — I,-- - ,4, corresponding to T(IS'), T(25), 
T(3S'), and T(4S'). We use the masses and widths 
from the Particle Data Group [13] ■ The widths 
are rT(is) 1.340(18) keV, rT(2S) = 0.612(11) keV, 
rx(3S) = 0.443(8) keV and T-t^as) ^ 0.272(29) keV. 
Given that the widths of the T(IS'), T(2S') and T(3S') 
were obtained at the same experimental facility, we 
will assume their uncertainties to be correlated. The 
masses are Mx(is) = 9.46030(26) GeV, Mr(2S} = 
10.02326(31) GeV, AfT(35) = 10.3552(5) GeV, and 
Mt(4S) = 10.5794(12) GeV. Finally, we use the effec- 
tive electromagnetic couplings from ^1^]. The BABAR 
Collaboration [l| has performed direct measurements of 
Rb in the continuum threshold region between 10.62 GeV 
and 11.21 GeV. There is also data on the full ratio R in 
the bottom-quark region by the CLEO Collaboration [l^ 
dating back to 1985. Subsequently, a later CLEO mea- 
surement in 1998 [i^, at a single energy, s ~ 10.53 GeV^, 
gives a total R-ratio roughly 30% lower than the 1985 
data in this region. Since this discrepancy remains unre- 
solved we shall use here only the BABAR data. As was 
pointed out in ^^l, these BABAR data cannot be used di- 
rectly in sum-rules, such as e.g. Eq.Q, for the following 
reasons. First, the initial-state radiation and the radia- 
tive tail of the T45 resonance must be removed. Second, 
the vacuum polarization contribution must be taken into 
account. We follow this procedure, as detailed in 0, to 
correct the BABAR data with results shown in Fig. [T] 
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FIG. 2: The values of mi,(10GeV), obtained for different 
values of sq and using the 10 different kernels in the class 
■Pg''"'"''^ (so, s). All results lie within the shaded region. 



tion C: The BABAR data are correct, and PQCD starts 
at y/s — 11.21 GeV. However, the PQCD prediction of 
Rhad is incorrect. The motivation for this option is that 
the exact analytical form of i?^^^^ is only known up to 
one- loop level. At order 0{a'^) already the full analytic 
result has to be reconstructed using Fade approximants 
to patch together information about n(s) obtained at 
y/s = 0, = 2TOh(/i) and y/s — >■ — oo. Both the Fade 
method, and the reliance on FQCD results obtained at 
threshold {^/s = 2mh{p)) could introduce unaccounted 
systematic errors. As a measure of the dependence of the 
method on the prediction of R^'^'~^^{s) up to sq (chosen 
to be large enough so that the high energy expansion 
becomes a rigorous prediction), we use <^^°(s) calcu- 
lated using the high energy expansion. The prediction 
of at = 11.21 GeV using the high energy ex- 

pansion is also closer to experiment than the prediction 
obtained using Rhad. 



The high-energy expansion of n(s), given in Eq.([5]), 
is only formally guaranteed to converge above y/s — 
4mf,(/i) « 15 GeV, due to non-planar diagrams having 
cuts starting there. Above this value the high energy 
expansion is an almost perfect approximation to the full 
analytic FQCD result [20]. Therefore, we shall always 
choose ^/so > 4mb(^) in Eq.([3]) so that it is safe to use 
the high energy expansion of n(s) in the contour integral. 
Between the end point of the data {^/s = 11.21 GeV) 
and ^Jso > iWibifi), we will use the best available FQCD 
prediction of Rb{s), obtained from the Fortran program 
Rhad [2^ . We consider this as data input, even though it 
stems from theory. The Rhad [23| prediction of Rb{s) is 
shown in Fig. [T] 

The first uncertainties affecting the bottom-quark 
mass are due to the uncertainty in the strong cou- 
pling as (Atts), the uncertainty in the experimental data 
(AEXF), and our limited knowledge of FQCD (A/i). The 
latter was estimated by varying the renormalization scale 

= 10 GeV by ± 5 GeV, running the mass calculated 
at this scale back to /i = 10 GeV and then taking the 
maximum difference. The second set are systematic un- 
certainties stemming from the fact that the FQCD pre- 
diction for Rb{s) does not agree with the experimen- 
tally determined values at the end point of the data 
{y/s — 11.21 GeV), as can be seen from Fig. [T] Two pos- 
sibilities for this discrepancy were considered in 01 . Op- 
tion A: The BABAR data are correct, but PQCD only 
starts at higher energies, say at ^/s = 13 GeV. Use then 
a linear interpolation between i?^-^^(l 1.21 GeV) = 0.32 
and i?^'^^°(13GeV) = 0.377, rather than the predic- 
tion from Rhad. Option B: The PQCD prediction from 
Rhad is correct, but the BABAR data are incorrect, per- 
haps affected by an unreported systematic error. In this 
case multiply all the data by a factor of 1.21 to make 
the data consistent with FQCD. In addition to these two 
options, we wish to consider a third possibility. Op- 



IV. CHOICE OF INTEGRATION KERNELS 



To minimize the dependence of results for the bottom- 
quark mass on Option A and Option C, the contri- 
bution from the region y/s = y/s* = 11.21 GeV to ^/sq 
should be quenched. This can be achieved by borrow- 
ing from the method of [2lj, where a Legendre polyno- 
mial was used to minimize the contribution of the then 
poorly known continuum threshold region. We choose 
here a Legendre-type Laurent polynomial, i.e. we con- 
sider linear combinations of powers of s chosen from the 
set S — {s^^, s^^, s~^, 1, s}. Inverse powers higher than 
lead to a deterioration of the convergence of PQCD, 
introducing large uncertainties from changes in the renor- 
malization scale fj, and the strong coupling as (see also 
[24|). We only use positive powers up to s^, as higher 
powers emphasize unknown 0{al) terms in the high en- 
ergy expansion. The optimal order of the Legendre-type 
Laurent polynomial was found to be 3 or 4. First, let us 
consider the order 3 case and let 

p{s) = V'i''-''^ (s, so) = A{s' + Bs' + Cs'=) , (10) 

subject to the global constraint 

/"So 

/ Vt''''\s,s^)s-^ ds^Q, (11) 

where n £ {0,1}, i,j,k G {—3,-2,-1,0,1}, and i,j,k 
are all different. The above constraint determines the 
constants B and C. The constant A is an arbitrary 
overall normalization which cancels out in the sum rule 
Eq.(I31). The reason for the presence of the integrand s~" 
above is that the behavior of Rb{s) in the region to be 
quenched resembles a monotonically decreasing logarith- 
mic function. Hence, an inverse power of s optimizes the 
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TABLE I: Results for mb(lOGeV) using kernels p{s) selected for producing the lowest uncertainty. Results from the kernels 
p{s) = and p{s) = used in 0|-|3| ^re given here for comparison. The errors are from experiment (AEXP.), the strong 
coupling (Aus) and variation of the renormalization scale by ±5GeV around fi = lOGeV (A/i). These sources were added in 
quadrature to give the total uncertainty (ATOTAL). The option uncertainties A A, AB and AC are the differences between 
m5(10 GeV) obtained with and without Option A, B, or C. As in these are not added to the total uncertainty, and are 

listed only for comparison purposes. 



quenching. As an example, taking sq = (16 GeV)^ (and 
A = 1) we find 

vt^-~^'°\s,so) = - (1.02 X 10-4 GeV-4) 

+ 3.70 X lO-'^ GeV-^ , (12) 

with s in units of GeV^. There are ten different kernels 
'p^,hk) ^ and the spread of values obtained for m;, using 
this set of different kernels will be used as a consistency 
check on the method. Outside the interval s S [s*,so]i 
•p(«j,fe)^g^ g^-j will blow-up, which leads to a suppression 
of the continuum threshold region relative to the well 
measured T-resonances. This will minimize the depen- 
dence of the results on Option B. Hence, this kernel min- 
imizes all three sources of systematic uncertainty. The 
fourth-order Laurent polynomial v''l'^''''^\s, s^) is also 
defined by the constraint Ea.([TT|). but with n £ {0, 1,2}. 

There arc also five different kernels V^''''^'^\s, sq). In 
general, the higher the order n oiVm the better the con- 
trol over the systematic errors. However, the price to 
pay is a reduction in the rate of convergence of PQCD, 
though this convergence can be improved by increasing 
Sq. In the Appendix we give explicit expressions for the 
various kernels used in table I. 



V. RESULTS AND CONCLUSIONS 



We considered a total of 15 different kernels p{s) used 

in Eq.([3]), 10 from the class of kernels 'P'^'^'^\s, sq) and 

5 from the class 'p4*'"'''^''^''(s, so)- All these are similarly 
constructed (i.e they obey Ea. ([TTt ). and hence have a 
similar ability to reduce the dependence of the bottom- 
quark mass on Options A, B, C. They do, however, 
place very different emphasis on theory. In particular, if 



say V^'''''^\s, So) only included inverse powers of s, then 
almost the entire right hand side of Eq.(|21) would em- 
anate from the residue, and hence from the low energy 
expansion of PQCD. If however P^*'^'''=^(s, sq) were com- 
posed of only positive powers of s, then only the high 
energy expansion of PQCD would enter the right hand 
side of Eq. ([3]) . Different kernels can therefore lead to sig- 
nificantly different dependencies on the renormalization 
scale /i. Our philosophy is to choose those kernels pro- 
ducing the lowest total uncertainty. The results from 
these are displayed in Table HI We also plot in Fig. 2 
the range of values for mf,(10GeV) obtained using all 
of the 10 kernels in the class 'P^''''^\s, sq), as a function 
of Sq. Remarkably, between 12 GeV < ^/so < 28 GeV, 
all of the masses obtained using all 10 kernels from 
the class V^~^'~^'^\s, sq) fall in the range 3621 MeV < 
TOfc(lOGeV) < 3625 MeV. Our method gives a consistent 
result even in the region .^so < 4mfc(/i) « 15 GeV where 
the high-energy expansion used in the contour integral in 
Eq.(|3]) is not guaranteed to converge. Using, rather, the 5 
kernels in the class v!^"'^''''^\s, sq), and varying sq in the 
range 18 GeV < ^/sq < 70 GeV, all of the masses thus 
obtained lie in the interval 3620 MeV < TOh(lOGeV) < 
3626 MeV. These results show a great insensitivity of our 
method on the parameter so, and also on which powers of 
s are used to construct v'^''''^\s, so) and V^'''''^'^\s, sq). 
This in turn demonstrates the consistency between the 
high and low energy expansions of PQCD. 
For our final result we choose the optimal kernel 
V^~^'~^''^\so, s) to obtain 

mf,(10 GeV) = 3623(9) MeV , (13) 

mb{rnb) = 4171(9) MeV . (14) 

This result is fully consistent with the latest lattice value 
TOh(lOGeV) = 3617(25) MeV [H. It is also consis- 
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tent with a previous QCD sum rule precision determi- 
nation giving TOf, (10 GeV) = 3610(16) MeV. Apart 
from our novel QCD sum rule approach, the inputs in 
the latter are almost identical to ours, with the excep- 
tion of their use of kernels of the form p{s) = s"", 
n G {2, 3, 4, 5}, and the use of a value of the strong cou- 
pling with a larger uncertainty. Their final result was 
obtained using p{s) = s~^, which can be seen from Ta- 
ble U as being far more sensitive to possible systematic 
uncertainties arising from Options A, B, C. They also 
determined mi, using p{s) — s~^, for which they obtained 
TOfc(lOGeV) = 3619(18) MeV. This value is closer to our 
result, which may not be surprising given that it is less 
sensitive to Options A, B, C thanp(s) — s~^, although 
not as insensitive as using our kernels. 
In conclusion, we have discussed here a finite energy QCD 
sum rule method with integration kernels involving in- 
verse and positive powers of the squared energy. The 
result for the bottom-quark mass has a lower total un- 
certainty, and is far less sensitive than the popular in- 
verse moment method to the three systematic uncertain- 
ties identified earlier, i.e. Options A, B, C. It should be 
appreciated from Table I that the results Eqs. (IT51) - P^ 
are independent of the PQCD prediction from Rhad in 
the region between ^/s ~ 11.21 GeV and ^/s — 4:rni,{^J,). 



VI. APPENDIX 



Up to an overall constant, the integration kernels 
'Pn{s,so) can be obtained from Eq.(ll). For complete- 
ness we list below the explicit expressions for all the poly- 
nomials used in Table I, at the corresponding values of 
So. First, for sq = (16GeV)2 

vt'^'~^'°\s,so) = - (1.015 X 10-^ GeV-^)s-i 
+ 3.694 X lO^^GeV"*^ , (15) 
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